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In this paper we derive an analytic approximation to the density of states for atoms in a combined optical 
lattice and harmonic trap potential as used in current experiments with quantum degenerate gases. We compare 
this analytic density of states to numerical solutions and demonstrate its validity regime. Our work explicitly 
considers the role of higher bands and when they are important in quantitative analysis of this system. Applying 
our density of states to a degenerate Fermi gas we consider how adiabatic loading from a harmonic trap into the 
combined harmonic-lattice potential affects the degeneracy temperature. Our results suggest that occupation of 
excited bands during loading should lead to more favourable conditions for realizing degenerate Fermi gases in 
optical lattices. 



I. INTRODUCTION 

Tremendous progress has been made in the preparation, 
control and manipulation of Fermi gases in the degenerate 
regime Q S S HBL S 0] ■ Such systems have many potential 
applications in the controlled study of fermionic superfluidity 
and the production of ultra-cold molecules. Another area of 
developing theoretical interest is in the physics of fermions 
in optical lattices (H |H [T(| [U d], and experiments have 
already begun to examine the properties of Fermi-gases (pre- 
pared as boson-fermion mixtures) in one-dimensional Il3lll4ll 
and three-dimensional UHfl^O optical lattices. 

Optical lattices have many features in common with crys- 
tals where a periodic lattice is also present, and many of the 
ideas and techniques from solid state physics have been ap- 
plied to this system. A unique property of optical lattices, as 
realized in experiments, is that the periodic lattice is accompa- 
nied by a harmonic confining potential, arising from an exter- 
nal magnetic trap or from effects related to the focused laser 
beams used to form the lattice. We shall refer to this potential 
as the combined harmonic-lattice potential, which is the main 
subject of the investigation presented in this paper (see Fig. 

LB. 

Whilst the harmonic potential and translationally invari- 
ant periodic potential are well characterized individually, their 
properties when combined are not as well-understood. Even 
though the harmonic trap is often much weaker than the con- 
finement provided by each lattice site, its effect on the spec- 
trum and properties of the system can hardly be considered 
small: it breaks the translational invariance of the system and 
changes the nature of the energy states in the deep lattice from 
compressed bands, to a set of unbounded overlapping bands. 
Several recent articles have considered aspects of this system 
d d d |H. In the context of a tightbinding model of 
ultra-cold bosons, the spectrum of the combined potential ap- 
pears to have been first considered by Polkovnikov et al. 12111 . 
Refs. B18lll9ll have made detailed studies of the combined po- 
tential spectrum (also within a tightbinding description), and 
closed-form solutions to this problem were recently given by 
Rey et al. [20]. In Refs. 11121 12211 an ideal gas of fermions 
in a ID combined potential was examined without making the 
tightbinding approximation. All of these studies have con- 



Combined Harmonic-Lattice 




Figure 1 : Schematic diagram comparing the translationally invariant 
lattice to the combined harmonic lattice considered in this paper. 



firmed that, for appropriate parameter regimes, parts of the 
single particle spectrum will contain localized states. This is 
in contrast to the translationally invariant system, where inter- 
atomic interactions or disorder are needed for localization to 
occur (e.g. see l23l l24ll ). In the combined potential local- 
ization arises solely from single-particle effects. Experiments 
with ultra-cold (though non-condensed) bosons 12511 have pro- 
vided evidence for these localized states. 

Many of the physical phenomena that are suitable to exper- 
imental investigation in optical lattices are sensitive to tem- 
perature and it is therefore of great interest to understand how 
the temperature of a quantum degenerate gas changes with 
lattice depth. Experimental results by Kastberg et al. j26ll 
in 1995 showed that loading laser cooled atoms into a three- 
dimensional optical lattice caused the atoms to increase their 
temperature 14111 . In previous work we have studied how 
the depth of a translationally invariant lattice affects the ther- 
modynamic properties of quantum degenerate Bose If27ll and 
Fermi [28] systems (also see ll29ll30ll3ll[32ll ). The most im- 
portant predictions of those studies relate to the temperature 
changes induced by the lattice, in particular that in appro- 
priate regimes increasing the lattice depth could be used to 
cool the system. For Fermi gases this cooling effect was used 
to predict that loading into an optical lattice could be used 
to enhance the conditions for observing the superfluid tran- 
sition IH]. However, recent work with a tightbinding model 
has shown that including the effects of the harmonic potential 
leads to an increase in the temperature of the system (com- 
pared to the Fermi temperature) by a factor of 2 13311 . demon- 
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strating the importance of fully treating the combined poten- 
tial. 

In this paper we consider the thermal properties of a Fermi 
gas in the combined harmonic-lattice potential including the 
role of higher vibrational states. We derive an analytic approx- 
imation to the single-particle density of states and compare 
this analytic density of states to numerical solutions and deter- 
mine validity criteria. We give characteristic energy and atom- 
number scales which can be used to predict when higher band 
effects will be important. Applying these density of states 
to a degenerate Fermi gas we consider how adiabatic loading 
from a harmonic trap into the combined harmonic-lattice po- 
tential affects the degeneracy temperature over a wide param- 
eter regime. Our results show that beyond tightbinding effects 
provide an appreciable correction to the adiabatic loading cal- 
culations presented in Ref. l33ll . Also, in the regime where 
higher band states are appreciably occupied, we find that the 
degeneracy temperature tends to increase to a lesser extent or 
even reduce during loading, suggesting that this regime might 
be useful for preparing a strongly degenerate Fermi gas in 
an optical lattice. Recently molecule production was used to 
measure the temperature of fermions in an optical lattice lfl5tl 
and using this (or other approaches B42I0 the predictions of 
our work can be examined in current experiments. 



II. FORMALISM 



difficult to arrange experimentally. However, our primary in- 
terest is in the thermodynamic properties of the system which 
are insensitive to this symmetry. 

We are interested in the limit where the lattice dominates 
the short-length scale properties of the system, and will take 
a <C aho, where ciho = min{ yj h/mujj }j =x ,y.z is the smallest 
harmonic oscillator length. This removes our need to consider 
systems where extremely tight harmonic confinement causes 
all the atoms to coalesce to a single site. 



1. ID spectrum 

Viverit et al. f\% have shown for the one-dimensional case 
of Eq. (fl} that when the lattice is sufficiently deep the eigen- 
states are localized to lattice sites. For our symmetric potential 
the eigenstate localized at site n will necessarily also localize 
at site — n, so strictly we should not call such states local- 
ized. However, this property is fragile to any asymmetry in 
the system, and has no discernible effect on the energy spec- 
trum (which is of primary interest to us). In this regime, the 
lattice site index n forms a convenient quantum number for 
the eigenstates, specifying the site where the state is local- 
ized. The respective energy eigenvalue is given by the value 
of the harmonic potential at that site, i.e. 
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A. Analytic density of states in the inhomogeneous lattice 
potential 

Here we consider the properties of a system described by 
the single particle Hamiltonian 



H 



2m 



V,. 



(1) 



where the combined potential V c is formed by an optical lat- 
tice potential and harmonic trap potential, i.e. 



Vc 



V h + Vu 

^m{u 2 x x 2 + uj' 2 y 2 +uj' 2 z z 2 ), 
V [sin 2 (fcir) + sin 2 (fcy) + sin 2 (fcz)]. 



(2) 
(3) 
(4) 



The harmonic trap is taken to be anisotropic, with angular fre- 
quencies {cu x , 0J y , 0J Z } along the coordinate directions. The 
lattice is of depth Vq, k is the wavevector of the counter- 
propagating light fields used to form the lattice and a = ir/k 
is the direct lattice vector. We also use k to define the re- 
coil frequency uir = Hk 2 /2m, and associated recoil energy 
Er = hjJR- The combined harmonic-lattice potential, as de- 
fined in Eqs. ©-(HJi has minima (i.e. lattice sites) at posi- 
tions r = n x ax + n y ay + n z az, where {x, y,z} are unit 
vectors, and {n x , n y , n z } are integers that are convenient for 
labelling particular lattice sites. We note that our highly sym- 
metric choice of the lattice potential, having a site coincident 
with the minimum of the harmonic potential, may be rather 



where lu is the trap frequency and we have set to zero the 
zero-point energy associated with the confinement in each lat- 
tice site. We will refer to these states as ground band states 
for clarity. The condition for localization is AE(n) > Jo, 
where AE(n) = e^°L — e^ ' is the difference in (harmonic 
trap) potential energy between lattice site n and n + 1, and Jo 
is the tunneling between sites [43]. This requirement is most 
difficult to satisfy near the trap potential minimum, where the 
difference in energy between adjacent sites is least. This va- 
lidity condition is equivalent to 



\n\ > n crit = q/8, 



(6) 



where n cl j t is the approximately the index of the closest site 
to the origin for which the localization condition is satisfied 
and q = 8Jo/mu> 2 a 2 is a dimensionless parameter which we 
discuss below (originally defined in 02OI1 ). 

While states satisfying Eq. (f5]) are valid for sufficiently 
large values of n, this expression neglects the existence of ex- 
cited vibrational states, which for the case of a translationally 
invariant lattice would correspond to the first excited band. 
The energy scale for the emergence of these excitations is 
e gap , and an analytic approximation for this quantity is cal- 
culated in Appendix [A] Like the ground band states, these 
states will also localize when the difference between potential 
energy at neighbouring sites exceeds the tunneling matrix el- 
ement for the first excited band, which we denote J\ . Where 
this condition is satisfied the spectrum of these states will take 
the analytic form 
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For clarity we shall refer to these as the first excited band 
states. This argument could be extended to additional bands 
of states, in particular, at an energy of roughly 2e gap the next 
vibrational states will be accessible. However, the tunneling 
rate Jm increases with the band index to, making the localiza- 
tion condition more difficult to satisfy, and analytic estimates 
for the energy gap to higher bands less accurate. Indeed, for 
single particle energies large compared to the lattice depth the 
spectrum will cross over to that of a harmonic oscillator. Ad- 
ditionally, for typical experimental parameters, the first two 
bands include sufficiently many states, and a large enough en- 
ergy range to provide an accurate description description of 
the system. Of course, high temperatures, or large lattice con- 
stants would require consideration of additional bands. 



2. 3D spectrum and density of states 
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For more than one spatial dimension the wavefunction lo- 
calization may be broken by neighbouring sites which have 
approximately the same local potential value. However since 
Eq. (fTJ is separable, the 3D ground band spectrum is com- 
pletely determined by the one dimensional results (0 and (0, 
and under the assumption of localized states (see below), the 
ground and first excited band spectra are given by the expres- 
sions 



e (0) 

,(1) 
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-ma (uj x n x + Lu y n y + u; z n z ) + e gap , (9) 



Figure 2: Energy level spacing for the ID eigenspectrum relevant to 
the experimental setup in Ref. 11611 . (a) Vo = 5Er, lo = 0.024lj_h 
(b) Vo = WE R , u = 0.033wr. Parameters derived for 40 K with 
lattice made from counter propagating A = 826nm lasers with har- 
monic confinement arising from the focused lasers (beam waist taken 
to be 50/ira). Labels A-E explained in text. 



discussion of the validity conditions of the spectra, additional 
bands will become accessible at e ~ 2e gap , and so Eq. ( TTOb 
should only be used in situations where ksT and the Fermi 
energy are much less than 2e gap . 



respectively. Because of the lattice symmetry, there are 3 
equivalent (and completely overlapping) first excited bands. 
Thus each energy state specified by the quantum numbers 
{n X) n y , n z } in (0 is three-fold degenerate [44] . This degen- 
eracy would be broken if the lattice depth was different in each 
direction causing the first excited bands to separate in energy, 
but we do not consider that case here. The separability of the 
potential means that the validity conditions discussed for the 
one-dimensional case apply immediately, using the respective 

(i) 

trap frequency in each direction, i.e. rij > n crit j — x, y, z, 

where n^L is n cr ;t evaluated according to (|6]l using the trap 
frequency along direction- j. 

The density of states for the spectra in Eqs. © and (O is 
given by 



9c(e) 
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where 6{x) is the unit step function, and ui = &uj x uj y ui z . The 
first term, corresponding to the ground band contribution to 
the density of states, exhibits a -yi-scaling with energy, sim- 
ilar to that of a homogeneous gas of free particles 113311 . The 
second term includes the contribution of the first excited band 
states that occur at energies e > e gap . As was noted in our 



B. Numerical Results for single particle spectrum 

1. ID spectrum 



Here we show some typical results of the ID spectrum of 
the combined potential in Fig. El In Refs. ifTsL IT9I [20I a 
detailed analysis of the spectrum has also been made, but in 
the tight binding (Hubbard) limit where only the vibrational 
ground state of each lattice site are included. Those studies 
considered a wide parameter regime of trapping frequencies 
and lattice depths, however in the tight-binding limit a sin- 
gle parameter describing the ratio of tunneling to harmonic 
potential is sufficient to characterize the nature of the eigen- 
spectrum. Several choices of parameter are used in the liter- 
ature, and we follow the choice of Rey et al. 12011 who de- 
fine q = AJq/ (^mu} 2 a 2 ). With the inclusion of higher bands 
this single parameter by itself is insufficient to characterize the 
spectrum and both the lattice depth and harmonic confinement 
parameters are independently important. In Fig. [2] we show 
the level spacing obtained from numerical diagonalization of 
the one-dimensional case of Eq. (Q~|). The parameters for this 
calculation were taken to correspond to those of the experi- 
ment |Q]1. For Fig. Ha) with V a = 5E R we find that q w 193 
and for Fig. Hb) with V = 10E R we find that q w 28. It 
is useful to compare our results to Fig. 1(b) of Ref. IU9I1 to 
assess the effects of higher band states on the spectrum. We 
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Figure 3: The critical site number for localization, as defined in Eq. 
{6j, as a function of lattice depth and harmonic trap frequency. Con- 
tour lines are shown for n cr it = 1, 3, 5 and 7. 

have segmented our spectrum with vertical dashed lines and 
indicated the characteristic regions by the letters A-E, which 
we explain: In region A tunneling dominates over the offset 
between lattice sites and the eigenstates are delocalized (for 
these states Eq. (0 will be a poor approximation). In region 
B localized states emerge which are approximately degener- 
ate in energy (the spacing between every second eigenvalue is 
approximately zero). Because of the localizations (and appli- 
cability of Eq. (O), the energy spacing between degenerate 
pairs increases linearly throughout this region. Regions A and 
B, as identified in Fig. [2] are equivalent to those in Fig. 1(b) 
of Ref. 11911 . and in general we find that the tightbinding de- 
scription is in good qualitative agreement with our more gen- 
eral analysis. For the next regions (i.e. C, D etc.) the role 
of higher bands is essential and cannot be described within a 
simple tightbinding model. Region C has a similar appearance 
to region A, and consists of delocalized excited band states. 
The scatter in level spacing seen in C occurs because ground 
band states are also available in this energy range. In region 
D the excited bands state have become localized (similar to 
what happened in region B to the ground band states). In re- 
gion E we have reached a sufficiently high energy scale that a 
3rd band of states have begun to contribute. 

The harmonic confinement used in Fig. [2] originates from 
the dipole conferment provided by the focused lasers used to 
make the lattice[45], and as the lattice depth increases so does 
the strength of harmonic confinement. We note that for Fig. 
1(b) of [19] q w 13 x 10 6 which is many orders of magni- 
tude away from that found in current experiments l46ll . The 
number of states in region A roughly scales as ^fq [ 20] so that 
in experimentally realized lattices region A is quite small and 
the majority of ground band states are well-localized. 



2. Critical site index 

In Fig. [3] we present n cr j t © for a wide range of lattice 
depths and trap frequencies (note that we have used band 
structure calculations in a translationally invariant lattice to 
determine Jo for each value of Vq). We observe that n cr it is 



large for small trap frequencies and shallow lattices so that 
only neighboring lattice sites quite far from the harmonic trap 
minimum have a sufficiently large potential difference to tun- 
neling ratio to cause eigenstate localization. In such cases the 
analytic approximation for the eigenspectra given in Eqs. (0 
and (O, and density of states given in Eq. ( fTOb will not be 
valid, and the result of the full numerical diagonalization will 
be necessary. 

With increasing lattice depth the ground band tunneling ma- 
trix element decreases and the potential difference to tunnel- 
ing ratio increases. Thus we find that n cr ;t decreases with 
increasing Vq. Similarly, increasing the harmonic trap fre- 
quency also leads to a decrease in n cr i t . The particular value 
of n cr j t that justifies the use of our analytic density of states 
([Tol l depends on the parameters of system under consideration. 
If the system extends over N s lattice sites in each direction, 
then n cr ; t <C N s will be sufficient to ensure that the major- 
ity of the occupied states are well described by the localized 
spectrum. 

3. Density of states 

Equation ( fTOb for the density of states in the combined po- 
tential is one of the central results of this paper. In this section 
we present numerical results to confirm the validity regime of 
this expression. To do this we diagonalize Eq. (UJ) to obtain the 
single particle eigen-spectrum {e^} for various trap frequen- 
cies and lattice depths. For the purposes of comparison to the 
analytic results, it is useful to construct a smoothed density of 
states, defined as 

z ^ e Jt-Ae j 

giving the average number of eigenstates with energy lying 
within Ae of e. 

In Fig. [4]we compare the numerically calculated smoothed 
density of states against the analytic result g c {e) for various 
lattice depths in an isotropic trap of frequency uo = Q.QhluoR. 
Agreement between the analytic and numerical calculations 
is seen to improve as the lattice depth increases. We also ob- 
serve that at energies greater than approximately twice the gap 
energy the analytic and numerical results begin to differ more 
significantly as the contribution of additional bands become 
important (the gap energy in each case is the energy at which 
the cusp in the analytic density of states occurs). 

It is of interest to more closely examine the reliability of the 
analytic density of states at low energy scales and for weak 
harmonic traps. In Fig. [5]we show such a comparison for an 
isotropic trap of frequency Co = 0.01^. For Vq = AEr the 
analytic result is in poor agreement with the numerical result, 
as the low energy states in this case are not localized. For 
Vq = 8Eh and more so Vq = 12.Br, the analytic result is 
seen to provide a useful description of the numerical density 
of states. To quantify the degree of agreement seen in these 
results it is useful to consider the energy of the (first) localized 
state at site n n cr j t , i.e. e cr it = m w 2 a 2 n 2 rit /2, indicated as 
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Figure 4: Density of states for uu = 0.051lor. Numerical results 
for the density of states are shown using various markers (as labeled 
in the figure). The analytic density of states given in Eq. \\0\ is 
plotted for the cases with V > 5-E.r (solid line), using the expression 
e gap = 2t/Vo Er — Er for the band gap energy. The harmonic 
oscillator density of states is also shown (dashed line). The averaging 
interval used for the smoothed density of states is Ae = 0.303-Er. 



and hence the lowest energy at which spectrum © is valid. 
For the case of Vq = 4E^, e cr it exceeds e gap so that no states 
in the ground band are localized below the energy scale at 
which excited band states become accessible. For the case 
of Vq — 8Er, e cr it — 0-5Er and for energies above this 
better agreement between the analytic and numerical results 
is observed in Fig. [5] For Vq = 12 En the critical energy is 
e cr it ~ 0.07-Er and the analytic density of states furnishes 
better agreement at much lower energies. 

In summary, we make the following observations about our 
analytic expression for the density of states: 

1. Within its regime of validity, the analytic density of 
states provides an accurate description at intermediate 
energy scales, i.e. at energies above e cr it where the 
spectrum is well-localized, yet below e ~ 2e gap where 
additional bands become accessible. 

2. The agreement between the analytic density of states 
and numerical calculations improves with increasing 
lattice depth and increasing trap frequency. 

III. APPLICATIONS TO FERMI GASES 
A. General Properties 

1. Fermi energy 
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Figure 5: The density of states for a combined harmonic trap and 
lattice in the loose trap limit. Isotropic harmonic trap with trap 
frequency u = O.OIujr. The critical energy, given by e cr it = 
mw 2 a 2 n 2 rit /2, is shown as a dashed vertical line and labeled by 
its corresponding lattice depth. The averaging interval used for the 
smoothed density of states is Ae = 0.10i?R. 



vertical dashed lines in Fig. [5] The quantity e cr it is the energy 
scale above which the ground band states become localized, 



The Fermi energy in the combined harmonic-lattice poten- 
tial, cf,c> 1S determined from the number of particles in the 
system according to 



N 



deg c (e). 



(12) 



To compute the Fermi energy from our analytic result ( fTOb it is 
convenient to define the cusp number l47ll N c as the number 



of particles for which 6f.< 



Because the energy gap 



depends on the lattice depth, so does N c 
( TT21 ) we obtain 



Using Eqs. ( fTOb and 



32 (ujr 
3tt 2 V w 
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(13) 



For N < N c (i.e. e < e gap ) only the first term in the density of 
states ( fTOb is nonzero in Eq. (fT2l) and we can invert to obtain 
the Fermi energy 



/ 3t: 2 N 
V 32 



2/3 



hw 2 



N <N C 



(14) 



For N > N c excited band states contribute. In this case a 
general analytic expression for the Fermi energy in terms of 
iV is not available, but for (N — N c ) -C iV c we make a series 
expansion of the integral of the density of states about e gap . 
This expression can then be solved perturbatively to yield the 
approximate expression 



egap 



1 + 5- 



2S 3 / 2 



1 + 3S 1 / 2 



0<<5<1 (15) 
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where S = 2(N — N C )/3N C is the small parameter. 

It is also convenient to define the cusp depth V c , as the lat- 
tice depth at which e FtC — £gap- The cusp depth is a function 
of the number of particles and harmonic trap frequency, given 
by 



Vc = 



hu>i 



3tt 2 N 
32 



2/3 



— +1 



(16) 



where we have made use of the analytic expression relating 
e gap to Vb, as derived in Appendix lAl 

The two cusp parameters characterize the interesting fea- 
tures of the combined harmonic-lattice system and can be in- 
terpreted as follows: 

• N c : For a system with fixed combined potential (i.e. 
fixed Vo, {ojj}), N c is the maximum number of atoms 
that can be accommodated in the ground band only. For 
N > N c the T — ground state of the system will 
contain excited band states. 

• V c : For a system with fixed atom number and harmonic 
trap frequencies, V c is the smallest lattice depth for 
which the atoms can be accommodated in the ground 
band. For Vq < V c the T = ground state of the sys- 
tem will contain excited band states. 

These parameters motivate us to emphasize the distinctive 
properties of the energy spectrum in the combined potential 
as compared to the usual periodic lattice case. In the transla- 
tionally invariant lattice, there is a fixed number of single par- 
ticle states in each band (equal to the number of lattice sites), 
and for sufficiently deep lattices (typically Vq > 2Er) the 
ground and first excited bands occupy disjoint energy regions 
separated by a finite energy gap. In contrast, for the combined 
harmonic-lattice potential, the energy bands are overlapping 
and can only be differentiated by the local nodal structure of 
the wavefunctions at each site, where they have a spatial char- 
acter approximately given by harmonic oscillator states (see 
Appendix IA1. This local structure of the wavefunctions is ap- 
parent in experiment, and leads to states of different bands 
(as we have defined them here) residing in distinctive regions 
of momentum space in expansion images (e.g. see Refs. 
lfl6i l36Tl ). If it is desirable to restrict the system to access 
only states of the ground band, so as to realize a system well- 
described by a Hubbard model, then according to our above 
prescription this necessarily requires N < N c or equivalently 
Vq > V c , in addition to having sufficiently low temperature. 



B. Sommerfeld analysis of isentropic loading of a degenerate 
Fermi gas 

The properties of a quantum degenerate Fermi gas can 
be approximated by the Sommerfeld expansion (e.g. see 
ll37tl ). Of particular interest is the expression for entropy 



S = ^-g(e F )k%T, valid for T < T F , where g(e F ) is the 
density of states evaluated at the Fermi energy, e F . For many 
applications to Fermi gas experiments, the parameter of most 
interest is the degeneracy parameter t = ksT/e F , i.e. the ra- 
tio of the temperature to the Fermi temperature, for which the 
Sommerfeld expression can be written as 



S = —e F g{e F )k B t. 



(17) 



Here we consider the change in degeneracy temperature of 
a Fermi gas as it is slowly loaded from a harmonic trap into 
the combined harmonic-lattice potential, as is done in exper- 
iments. To characterize this temperature change we assume 
that the loading is isentropic so that the initial entropy in the 
harmonic potential (with initial degeneracy temperature U) is 
the same as the final entropy when the system is in the com- 
bined harmonic-lattice potential (with final degeneracy tem- 
perature tf). Within the validity regime of the Sommerfeld 
relation ( fTTj l, the ratio of these temperatures is given by 



U 



K C (N) ' 



(U,t f < 1) 



(18) 



obtained by assuming that S remains constant, where we have 
introduced the dimensionless extensive parameter k x (N) = 
e F , x g x {e Fx ) (x = h,c), with g h (e) and e Fth the density of 
states and Fermi energy for a harmonic trap respectively. We 
have chosen to express k as a function of N rather than e F , 
since the number of atoms remains constant during the load- 
ing procedure, whereas the Fermi energy may change signifi- 
cantly. 



For the purely harmonic trap gh{e) 



e F<h = frw(6A0 1/3 ( e -g- see H), which give 

K h {N) = 3N. 



V(2ft 3 



and 



(19) 



For N < N c , the excited band states do not contribute to k c 
and using results ( TTOb and ( TT4l > we obtain 



K C (N) 



-N, 



N <N C . 



(20) 



Thus for N < N c the ratio of degeneracy temperatures will 
increase by a factor of two. Because 7V C is a monotonically 
increasing function of Vb, it might be expected that for suf- 
ficiently deep final lattice depth we will always obtain this 
factor of 2 increase in the degeneracy temperature. However, 
N c decreases with increasing ui, and since Q may change with 
lattice depth (e.g. if the lattice is produced by focused lasers) 
a large final lattice depth may in fact lead to N > N c . We 
also note that when N < N c the ratio of the degeneracy tem- 
peratures ( fT8l is independent of the harmonic trap frequency, 
even if this parameter changes during the loading [48]. Simi- 
lar conclusions to those presented in this section are given in 
Ref. JH. 
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Figure 6: Temperature of an ideal Fermi gas isentropically loaded into a combined harmonic-lattice potential. The ratio of the final degeneracy 
temperature (after loading) to the initial degeneracy temperature in the harmonic trap are shown for various initial temperatures and lattice 
depths for the case (a) TV = 50 x 10 3 and (c) TV = 250 x 10 3 . In (e) the Fermi energy is shown as a function of lattice depth for TV = 25 x 10 
(circles), TV = 250 x 10 3 (crosses). (b), (d) and (f) correspond to (a), (c) and (e) respectively, but are calculated using the analytic density of 
states J 10b . For all results we have used an isotropic harmonic trap of frequency of Q = 0.04ojr. 

C. Numerical results for isentropic loading of the combined harmonic-lattice 

I 

In this section we numerically examine the temperature of potential for a range of lattice depths. Our main results are cal- 
an ideal Fermi gas loaded into a combined harmonic-lattice culated using the energy spectrum found by numerically diag- 
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onalizing ([T). To determine the temperature under this type of 
loading we calculate the entropy of a gas of N fermions over a 
range of temperatures and lattice depths, i.e. S(T, Vo, C3, N), 
where uj — {u) x ,u> y ,u z } . We numercially invert this func- 
tion to find temperature as a function of the other quantities 
T(S, Vq,uj, N), and by examining the behaviour of T for fixed 
S, we can predict the temperature of the gas as a function of 
lattice depth. 

Our procedure for determining entropy is as follows: The 
single particle spectrum {ej} of the lattice is calculated for 
given values of Q and Vo- We then calculate the partition func- 
tion Z 

logZ = ^log(l + e-' 3 ( £ J-^) , (21) 

3 

where fi is found by ensuring particle conservation. The en- 
tropy of the system can then be expressed as 

S = k B (log Z + (3E - npN) , (22) 

where (3 = 1/ksT, and E = —d\nZ /d(3 is the mean energy. 

In Figs. [6ja)-(f) we show the properties of an isentropically 
loaded gas for various parameters. In Figs. [6ja) and (c) we 
show the ratio of the final to initial reduced temperatures for 
N = 50 x 10 3 and N = 250 x 10 3 respectively. In both 
cases, the reduced temperature is seen to increase as the lat- 
tice depth increases. For N = 50 x 10 3 atoms the critical 
lattice depth is V c « 2.3Er. For the lowest initial reduced 
temperature (i.e. the front-most curve in|6ja)) we see that the 
reduced temperature increases by a factor of two by the time 
that Vq increases beyond V c , in agreement with the predic- 
tions of Eqs. (fT8l-d20t (also see Ref. JH). Similarly, for 
the case of N — 250 x 10 3 atoms, V c w \2AE R and for 
the lowest temperature result in Fig. |6jc), we see that the re- 
duced temperature increases by a factor of two as Vq increases 
beyond this value of V c . For higher initial temperatures the 
Sommerfeld result does not hold. In Figs. |6ja) and (c) we see 
that the warmer systems (larger t{ values) have the contrast- 
ing behaviour of heating up more or less than the Sommerfeld 
prediction respectively. In that regime the degeneracy tem- 
perature is dominated by the change in the Fermi energy that 
occurs during the lattice loading procedure, as shown in Fig. 
|6je). That is, for the case in Fig. |6ja) |6jc)] the Fermi energy 
tends to decrease [increase] with increasing lattice depth. 

In Fig. |6jc) (and to a lesser extent in|6ja)) we see that while 
excited bands are occupied (i.e. for V < V c ) the degeneracy 
temperature increases by a factor of less than 2. This sug- 
gests that having excited bands occupied might provide more 
favourable conditions for investigating fermionic superfluidity 
in lattices. Additionally, because the tunneling rate is larger 
for higher bands it may be more easy to reversibly manipulate 
the lattice in this regime. 

In Figs. |6fb), (d), (f) we show the results equivalent to those 
in Figs. Ufa), (c), (e), but calculated using the analytic den- 
sity of states given in Eq. ( TTOb . Qualitatively the agreement 
between the results is good for Vq > 4Er. The main dis- 
crepancy is observed for small Vq values where the role of 
non-localized states and higher bands is important. 




Figure 7: N c for the case of 40 K in a lattice formed by focused lasers 
(solid) A = 826nm with 70[im waist, (dashed) A = 1200nm with 
50/im waist. 



IV. RELATION TO EXPERIMENTS 

It is of interest to compare how important excited band ef- 
fects might be for current experiments. In Fig. [7] we show 
N c for parameters similar to those used in recent experiments 
by the ETH group lfl5l [Trjll . In those experiments up to 10 5 
atoms were prepared in each spin state. Because the harmonic 
confinement increases with lattice depth (uj ~ y/VoEu/h), we 
find that the Fermi energy will eventually lie in the first excited 
band, but only for very large lattice depths (Vo > 280Er). 
We also consider a longer wavelength lattice made from lasers 
with A = 1200nm with slightly tighter focus (beam waist of 
50/im) for which the recoil energy and the gap to higher bands 
is smaller. In such a configuration (dashed curve) we see that 
higher bands would become important for much lower atom 
numbers. 

In Fig. [8] we consider the effect of loading on the degener- 
acy temperature for parameters relevant to the experiment in 
Ref. II 1511 . The results in Fig. [Sfa) are for the same parameters 
considered by Kohl who made a calculation in the tightbind- 
ing limit (see Fig. 1 of Ref. ll33ll ). We broadly find agreement 
with those results, however make note of several differences. 
First, at the lowest depths considered (Vq ~ 5 En) Kohl ob- 
served the reduced temperature to initially decrease with in- 
creasing lattice depth. In contrast our results, which are valid 
for all lattice depths, do not exhibit this feature, and instead 
are seen to smoothly connect with the Vq = QEr (purely har- 
monic case). We conclude that in this region the tightbinding 
approximation and the neglect of higher bands is not valid. 

Second, in Kohl's results the reduced temperature is ob- 
served to saturate to a value of twice that of the initial har- 
monic trap when the lattice is sufficiently deep (typically 
Vq > IHEr), Many of our high temperature curves instead 
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Figure 8: Degeneracy temperature for isentropic loading of 40 K into 
an optical lattice, (a) N = 50, 000 and (b) N = 750, 000. Parame- 
ters correspond to solid line in|7] Additional harmonic confinerment 
of frequency uj = 0.005o;jj is superimposed on the confinement due 
to the focused beam waist to make the spectrum well-behaved as 

Vb->o. 



show a slight decrease in reduced temperature as the depth 
of the lattice increases. We have verified that this behaviour 
is due to higher band states, and that if we neglect them from 
our calculations our results saturate in agreement with those in 
ll33ll . This indicates that even when the deep lattice behaviour 
is dominated by the ground band (i.e. we have N < N c ), if 
a small number of atoms are able to thermally access excited 
band states they can have a significant effect on the tempera- 
ture of the system during loading. 

In Fig. [Sib) the same parameters are used, except the num- 
ber of atoms is increased to 750, 000. According to Fig. [7] 
for this number of atoms V c « 15-Er. In this case we see 
that only a small increase in the degeneracy temperature oc- 
curs during loading and ultimately for sufficiently large final 



lattice depth (i.e. Vq > 2QEr) the reduced temperature is ob- 
served to be approximately the same as in the initial harmonic 
trap. We note that this result requires the occupation of higher 
bands and cannot be analysed using a tightbinding approach. 

The discussion in this section shows that while current ex- 
periments are likely not strongly affected by higher bands, 
the parameter regime where they become important is rather 
close. Experiments could enter this regime by using larger 
wavelength lattices or tighter harmonic confinement (e.g. 
more tightly focused lasers to produce the optical lattice). 



V. CONCLUSIONS 

We have discussed the nature of the spectrum of a Fermi gas 
in a combined harmonic trap and optical lattice potential. Us- 
ing this spectrum we have derived an analytic density of states 
that is relatively accurate for the lattice depths and harmonic 
confinements used in experiments. We have characterized the 
validity criteria for this density of states and have used it to 
characterise a Fermi gas in the combined potential. As an 
application we have examined how adiabatic loading from a 
harmonic trap into the combined harmonic -lattice potential af- 
fects the degeneracy temperature of an ideal Fermi gas. Our 
results show that when excited band states are occupied the 
system is less heated by the lattice loading, and may be less 
sensitive to non-adiabatic effects, suggesting that this regime 
is worthy of further investigation in experiments. 
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Appendix A: ANALYTIC APPROXIMATION TO BAND GAP 
ENERGY 

In this section we derive an analytic expression for the band 
gap in deep lattices. We go beyond the usual harmonic oscilla- 
tor approximation and obtain results equivalent to those used 
inRef. (H. 

We will consider the standard harmonic approximation the 
optical lattice potential, by making a Taylor expansion about 
the lattice site minimum at x = 0, i.e. 

V^^H-cosM]^^-^, (Al) 

where b — 2k is the reciprocal lattice vector. 

Casting the harmonic term in the form of a harmonic oscil- 
lator potential, ^moj^^x 2 , yields the effective harmonic os- 
cillator frequency of WLatt = \/Vo& 2 /2m, and the localized 
states in the optical lattice can be approximated as harmonic 
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Figure 9: Comparison between numerically determined energy gap 
for a ID sinusoidal lattice (solid), with the upper and lower limits 
corresponding to the band gap at the center and the edge of the Bril- 
louin Zone. The harmonic approximation, huihaxt = 1\JVqEr (cir- 
cles) and the estimate treating the quartic expansion perturbatively , 
e g ap = 2VVo-Er - E R (crosses). 



oscillator states. This approximation neglects the influence of 
tunneling between sites that gives a quasimomentum depen- 
dence, and predicts that the band gap is equal to /iWLatt- For 
our purposes it is desirable to go beyond this approximation 
and obtain a more accurate analytic expression for the band 
gap energy. To do this we use the harmonic oscillator states 



to treat the quartic term in Eq. ( IAU perturbatively. This is 
most easily done using the normal ladder operators, so that 

2m % - ) + a) 4 , which gives the 1st order shifts 
in the oscillator state energies as 



AE n = 



^(2n 2 + 2n+l). (A2) 



Thus the approximate energies of the localized states are 
given by 



1/2)- i(2n 2 + 2n- 



E R . (A3) 



Of most interest is the difference in energy between the n = 
and n = 1 states, which provides an estimate of the energy 
gap between the ground and first excited bands, i.e. 



e gap — 2y/VoEn — Ej 



(A4) 



We see that treating the quartic term leads to a 1 recoil sup- 
pression of the band gap compared to the harmonic oscillator 
frequency. In Fig. [9] we compare this analytical expression 
with the band gap determined numerically by evaluating the 
full band structure of the ^ [1 — cos(fox)] potential. For shal- 
low lattices the band gap is strongly dependent on the value of 
quasimomenta considered and error bars indicate the range of 
energy gaps from band centre to band edge. 
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